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ABSTRACT 

We report the results of the first two 5 GHz e-VLBI observations of the X-ray 
binary Cygnus X-3 using the European VLBI Network. Two successful observing ses- 
sions were held, on 2006 April 20, when the system was in a quasi-quiescent state 
several weeks after a major flare, and on 2006 May 18, a few days after another flare. 
At the first epoch we detected faint emission probably associated with a fading jet, 
spatially separated from the X-ray binary. The second epoch in contrast reveals a 
bright, curved, relativistic jet more than 40 milliarcseconds in extent. In the first, and 
probably also second epochs, the X-ray binary core is not detected, which may indicate 
a temporary suppression of jet production as seen in some black hole X-ray binaries in 
certain X-ray states. Spatially resolved polarisation maps at the second epoch provide 
evidence of interaction between the ejecta and the surrounding medium. These results 
clearly demonstrate the importance of rapid analysis of long-baseline observations of 
transients, such as facilitated by e-VLBI. 

Key words: accretion, accretion discs - stars: individual: Cygnus X-3 - ISM: jets 
and outflows - radiation mechanisms: non-thermal - techniques: interferometric. 



1 INTRODUCTION 

Th e X-ray binary Cygnu s X-3 was first dete cted in X-rays 
by Giacconi et all (|l967l ). T he infrared (e.g. Becklin et al.l 
and X-ray fluxes (e.g. iParsignault et al.lll972l ) show 
a periodicity of 4.8 hours which is interpreted as the 
orbital period of the sy stem. The nature of the com- 
pact object is not known (|Schmutz. Geballe fc Schiidll 19961 : 
I Mitral Il998f ). As for the companion star, there is com- 
pelling evidence po i nting toward a WN Wolf-Raye t star 
Jvan K crkwiik ct al. 1996; Fe nder. Hanson fc Poolevl 1 19991 : 
iKoch-Miramond et al.ll2002l ). 

Giant outbursts and large flares have be en observed at 
radio wavelengths in Cygnus X-3 since 1972 l)Gregorv et al.l 
1972). In quiescence the soft X-ray emission is correlated 
with the radio emission, while the hard X-ray is anti- 
correlated with the radio; in a flare state, the situation is 
reversed: the hard X-ray correlates with the radio and the 
soft X-ray emission is anti-correlated dWatanabe et al ] | 19941 ; 
iMcCollough et al J 1 19991 ; IChoudhurv et alj|2002h ~ 
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Radio observations made during such large flares 
at different resolutions with the Very Large Array 
(VLA), Multi-Element Radio Linked Interferometer Net- 
work (MERLIN), Very Long Baseline Array (VLB A) 
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or are consistent with two-sided relativistic jets (with the 
notable exception of the VLBA observations of a flare in 
February 1997, wh e n the jet was apparently one-sided; 
iMioduszewski et all (120011 )). 



2 OBSERVATIONS 

One of the aims of e-VLBI is to enable mapping with long- 
baseline networks of radio telescopes in a manner which 
makes it possible to map transient phenomena, such as mi- 
croquasars, in near real-time. This will provide the ability 
to make informed decisions about the optimum observing 
strategy to employ (frequency of observations, array compo- 
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sition, calibrator selection, etc.) and the need for repeated 
mapping observations, as well as greatly simplifying the ob- 
serving procedure and improving its reliability. In the case 
of the EVN, the data are transported to the correlator at 
JIVE (Joint Institute for VLBI in Europe) in real-time using 
IP routed networks, connecting the radio telescopes through 
national research networks and the pan-European research 
network GEANT 2 via the Dutch national research network 
SURFnet to JIVE. At the time of our observations, a sus- 
tainable data rate of 128 Mbit s _1 could be achieved. This is 
expected to improve significantly in the near future. Descrip- 
tio n of some of the deve l opment work in this fi eld was given 
by iParslev et al.l i|2003l '); ISzomoru et al] l|2004 ) and reports 
of more recent development are in preparation (Szomoru et 
al.; Strong et al.). The data were transferred usin g Mark5A 
recor ders with 1 Gbit s _1 network interface cards |Whitnevl 
l2003h . 

We observed Cygnus X-3 at 5 GHz on three occasions, 
on 2006 March 16, April 20 and May 18 using the six anten- 
nas forming the current European e-VLBI network (e-EVN): 
Onsala 20 m, Torun, Jodrell Bank Mkll, Cambridge, Medic- 
ina and Westerbork (phased array). The observing run on 
March 16 was made in response to the first open call for e- 
EVN observations two weeks earlier. Unfortunately, due to 
technical problems with the off-line correlator software, the 
data from March 16 were not useable and have been omit- 
ted from the present study. In our second e-VLBI session, 
on April 20, we observed Cygnus X-3 for an effective inte- 
gration time of about 7 hours in a complex experiment con- 
sisting of interleaved observations of Cygnus X-3 and GRS 
1915+105 (see the accompanying paper - Rushton et al. 
(submitted) - for the latter target). Occasionally, delay and 
rate jumps were introduced in the data at the correlator 
Station Units, caused by a combination of limited memory 
buffer size and long correlation jobs. Although the effect was 
station based and thus the closure phases were unaffected, 
this may have resulted in erroneous phase transfer in a cou- 
ple of phase-referencing scans. The problem was eventually 
solved by self-calibrating the target data. During this run, 
the target was in a 'quiescent' state, with the flux density 
levels in the decay phase following the major radio flares in 
2006 March. At the beginning of May, Cygnus X-3 under- 
went another significant flare and we observed it during an 
active state via a target of opportunity proposal for around 
12 hours on May 18 (Fig. 1). During both observations the 
amplitude of the flux density variations was around 10 per- 
cent, except for some ~1 h long larger enhancements in the 
flux density registered at the beginning and end of the April 
run and at the beginning of the May run. Although this 
might have introduced artifacts in the radio images, we are 
confident that they would not be of significant importance. 

The data were calibrated in AIPS (e.g. iDiamoii d Hl995T l 
using standard procedures. For phase-referencing we used 
J2007+4029, a bright (~2 Jy correlated flux) calibrator lo- 
cated 4? 7 away from the target, employing a cycle time of 
8 min (5 min on the target, 3 min on the phase-calibrator). 
The solutions derived from J2007+4029 were smoothed, ex- 
trapolated and applied to Cygnus X-3. The a priori am- 
plitude calibration was improved using the unresolved cali- 
brator J2002+4725. Bandpass calibration was employed, us- 
ing J2007+4029 in epoch I (2006 April) and J2002+4725 in 
epoch II (2006 May). We used J2007+4029 as fringe finder 
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Figure 1. Top: RXTE/ ASM hardness ratio HR2=5-12 keV / 3 - 
5 keV. Middle: RXTE/ ASM 2-10 keV X-ray light curve. Bottom: 
15 GHz radio (Ryle Telescope) light curve. Vertical lines indicate 
the dates of our e-VLBI observations. 

in the former run and OQ208 in the latter. The data were 
self-calibrated once in phase with a 30-s solution interval and 
once in amplitude and phase simultaneously with a 30-min 
solution interval. The resulting radio maps are presented in 
Fig. 2. 

The instrumental polarisation leakage (the D-terms) 
was also determined for epoch II from observations of the 
unpolarised calibrator OQ208. For calibration of the elec- 
tric vector position angles (EVPAs) at VLBI scales we ana- 
lyzed separately the J2007+4029 Westerbork Synthesis Ra- 
dio Telescope (WSRT) data. We found that this calibrator 
has a fairly compact polarisation structure at VLBI scales 
enabling us to calibrate the EVPAs with an estimated error 
of 10 degrees. The distribution of the EVPAs (in the observer 
frame) and the fractional linear polarisation of Cygnus X-3 
are shown in Fig. 3. The level of Faraday rotation due to the 
galactic interstellar medium, affecting the overall distribu- 
tion of EVPAs, is unknown and therefore unaccounted for. 
For an arbitrary high rotation measure in the g alactic plane 
|Simard-Normardin. Kronberg fc Buttonlll98ll ) of, say, 200 
rad m -2 , the corresponding rotation of the EVPAs at 5 GHz 
would be up to ~40°. 



3 RESULTS AND DISCUSSION 

Fig. 2 presents the images from 2006 April 20 (epoch I) and 
May 18 (epoch II). Epoch I reveals a single faint component, 
possibly slightly resolved, while epoch II reveals a bright 
curved structure based around three bright knots. Also in- 
dicated on Fig. 2 is the location of the X-ray binary radio 
core as estimated by Miller- Jones et al. (2004; consistent 
also with Mioduszewski et al. 2001). 

In order to facilitate comparison between these two 
epochs, and more generally to allow us to extract as 
much information as possible from the observations we fit- 
ted the data in the uv-p\a.ne with elliptical Gaussians us- 
ing DIFMAP (|ShepherdTl997M . For epoch I we used one 
Gaussian while for epoch II we used as a starting model 
three Gaussians at the positions of the bright features as 
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Figure 2. Left: 5 GHz e-VLBI radio map of Cygnus X-3 on 2006 April 20 (epoch I). The contours are at -15, 15, 25, 35, 45, 55, 65, 75, 
85 and 95 times the rms noise of 0.03 mjy beam -1 . The beam size is 7.9 X 7.4 mas 2 at PA=-19?6. The colour code bar on top of the 
map is expressed in mjy beam -1 . Right: 5 GHz e-VLBI radio map of Cygnus X-3 on 2006 May 18 (epoch II). The contours are at -15, 
15, 25, 35, 45, 55, 65, 75, 85 and 95 times the rms noise of 0.2 mjy beam -1 . The beam size is 7.6 X 6.9 mas 2 at PA=88?3. The colour 
code bar on top of the map is expressed in mjy beam -1 . The positions of the centres of the elliptical Gaussians fitted to the data are 
represented by a cross for epoch I and stars for epoch II. The estimated location of the radio core from Miller- Jones et al. (2004), J 2000 
coordinates: RA 20 h 32 m 25f 77335, Dec. 40°57'27''9650, is indicated with a diamond. Note the change in scale between the two images. 



seen in the radio image before self-calibration. The re- 
sults are summarized in Table 1. The errors in the coor- 
dinates consist of the errors in the position of the calibrator 
J200 7+4029 (a a =0. 000180 s; ^=0.00240 arcsec; iFev et all 
2004) plus the error in the phase-referencing procedure 
itself (estimated as 0.3 mas; iPradel, Chariot fc Lestradei 
2006). We note that due to the poor itt>-coverage at short 
baselines the fitted models are not well constrained, par- 
ticularly the size of the Gaussians reported are proba- 
bly slightly underestimated. Scatter-broadening is an is- 
sue for Cygnus X- 3; its geometric mean size is 15 ± 
6 mas at 5 GHz dWilkinson, Naravan fc Spencer! Il994l; 
ISchalinski et~ai1l 19951 ) and might be correlated with the flux 
at short VLBA baselines (|Newell. Garrett~fc Spencer|[l998l) . 

We draw attention to the fact that the maps shown in 
Fig. 2 account for a small fraction of the total flux density 
of the source. The 5 GHz e-VLBI integrated flux density 
measured in the image plane for epoch I is 4.8 mjy. How- 
ever, in the uv-plane, the data were fitted well by a (larger) 
elliptical Gaussian with a flux density of 87 mjy. For epoch 
II, the total flux density in the uv-plane (determined by 
adding the flux densities contributed by the three elliptical 
Gaussians fitted to the data) is 909 mjy (much of the flux 
being observed on the shortest baseline), while in the im- 
age plane only 200 mjy is recovered. The explanation for 
these differences resides in the fact that because of the lack 
of «i>-coverage at short baselines the large scale emission 
could not be imaged properly. It is also worth mentioning 
that the simultaneous Westerbork integrated flux density in 
epoch II is 1.5 Jy, so at the e-VLBI scale we are missing al- 
most 40 per cent of the emission due to the unavailability of 
very short baselines. As an overall comment, the relatively 
sparse sampling at intermediate baselines has resulted in 
some considerable loss of flux and uncertainty as to the true 



strengths of individual knots. While we are confident that 
the morphology of the imaged structure is real, it is clear 
that there is a significant, more diffuse component present 
which we were unable to image. 

The single radio-emitting component from epoch I and 
knot C from epoch II are positionally coincident within un- 
certainties. This identification of the radio emitting regions 
in quiescence and outburst might intuitively lead to the con- 
clusion that the core was located here. However, the large 
positional discrepancy with the previously estimated core 
position, plus the high degree of linear polarisation (10 per 
cent, Fig. 3) in the region suggests that knot C is probably 
not the core (see further discussion below). Therefore, epoch 
I probably shows not the core but extended radio emission 
physically separate from the X-ray binary which may be as- 
sociated with the 2006 March flaring events (Fig. 1). This 
in turn implies that the core itself was below detection lev- 
els, with a 3cr upper limit of 0.4 mjy, about two orders of 
magnitude fainter than typical levels for the source, and 
four orders of magnitude fainter than the brightest flares. 
Such 'quenched' emission, lasting for periods of weeks to 
mon ths, was reported before other major flares of Cygnus 
X-3 (|Waltman et al.lll994[l995l ; lFender et al.lll997l ). 

Epoch II reveals a much brighter and more complex 
structure, with the appearance of a curved jet with three 
bright knots. Knot A is clearly much closer to the previ- 
ously estimated core location, but is still ~24 mas away. Mio- 
duszewski et al. (2001) and Miller- Jones et al. (2004) report 
the same core position to within 1 mas in data separated 
by 4.5 years (but note that this is rather uncertain given 
the fact that Mioduszewski et al. reported that their phase- 
referencing was not e ntirely successful). F or a distance to 
Cygnus X-3 of -9 kpc (|Predehl et al.ll200QD and an arbitrary 
speed in the plane of the sky of 100 km s~ , the upper limit 
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Figure 3. 5 GHz e-VLBI radio map of Cygnus X-3 on 2006 May 18 (epoch II) with the distribution of the EVPAs (left) and of the 
fractional linear polarisation (right) i.e. the ratio between the polarised flux and the total I flux, (Q 2 +U 2 ) 1//2 I — 1 with I, Q and U being 
the Stokes parameters. The code bar is expressed in percentage. In both images the Stokes I contours are at -15, 15, 25, 35, 45, 55, 65, 
75, 85 and 95 times the rms noise of 0.2 mjy beam -1 . 



Table 1. Results of the nt>-planc elliptical Gaussian fitting. The formal errors in coordinates (see text) are 0.000200 s in RA and 0.00270 
arcsec in Dec. 



Epoch 


Component 


RA (J2000) 


Dec. (J2000) 


Size 


PA 


Flux density 










(mas X mas) 


(dag) 


(mjy) 


epoch I - April 20, 2006 


single 


20 h 32 m 25? 771437 


+40°57'27'.'92800 


17.5 ± 0.3 x 8.8 ± 0.5 


150 ± 2 


87 ± 9 


epoch II - May 18, 2006 


A 


20 h 32 m 25? 771628 


+40°57'27'.'96566 


23.2 ± 0.7 x 10.6 ± 2.5 


69 ± 13 


60 ± 6 




B 


20 h 32 m 25? 771051 


+40°57'27'/95621 


38.4 ± 2.1 x 30.7 ± 8.1 


82 ± 22 


748 ± 17 




C 


20 h 32 m 25 ! : 771525 


+40°57'27'.'93000 


19.4 ± 1.1 x 12.2 ± 2.3 


55 ± 7 


101 ± 3 



to the positional shift due to the proper motion in the time 
interval between the observations of Miller- Jones et al. and 
ours is ~10.5 mas. Higher speeds are definitely possible (for 
instance GRS 19 15+105, at ~11 kpc has a proper m otion 
of ~5.8 mas vr' MDhawan. Mirabel fc RodriguezhoOOl ) ■ but 
given the present knowledge we cannot confidently identify 
the proper motion of the system as the cause of the shift. It 
seems likely therefore that even knot A may not correspond 
to the core, but simply to the nearest knot to it, indicating 
that at this epoch the core is also quenched. The orientation 
of the knots with respect to the core implies that all three 
are associated with the same side of the jet (likely to be 
approaching, although the complex environment of Cygnus 
X-3 may strongly influence appearances), and were ejected 
in the sequence C - B - A. Associating knot C with the 
large radio flare which occurred about a week earlier would 
imply proper motions in the jet of order 10 mas d" 1 , consis- 
tent with those reported in Miller- Jones et al. (2004). For a 
distance of ~9 kpc, this corresponds to a projected velocity 
of ~ 0.5c, and an intrinsic jet velocity at least as large. 

It is interesting to compare the 'quenching' of the core 
radio emission to that observed in the soft X-ray state o f 
black hole c a ndida te X-ray binaries (e.g. lFender et al1ll999t ). 
iGallo et al.l (|2003l ) have reported a correlation between X- 
ray and radio luminosities in black hole candidates in the 
hard X-ray state, below which the radio emission falls dra- 



matically in softer X-ray states. Assuming a distance of 9 
kpc, the core radio emission of Cygnus X-3 in epoch I is 
almost two orders of magnitude below the correlation, com- 
parable to the strongest limits on quenching found for other 
systems. Since neutron stars do not appear to show such 
strong quenching in soft X-ray states (jMigliari et al.l (|2004l ) 
- admittedly this is based on a very small sample) , this may 
be circumstantial evidence that Cygnus X-3 contains a black 
hole (although the extreme nature of the environment in this 
system might be responsible in some other way). However, 
the extreme scattering of X-rays in the system is always go- 
ing to make it hard to relate the X-ray behaviour to that 
of the other black hole candidates. It is noteworthy that in 
many cases 'residual' radio emission, such as that we believe 
we are imaging in epoch I, will limit the extent to which ob- 
servations on arcsec or coarser angular scales can measure 
the true extent or rapidity of such quenching. As a result, 
VLBI observations such as these are required to truly probe 
rapidly-quenched radio states in bright jet sources. 

In epoch II, the integrated flux densities from WSRT 
(1.5 Jy at 5 GHz) and the Ryle Telescope (0.8 Jy at 15 GHz; 
Fig. 1) reveal a steep spectrum (oc v~ 0,6 ), indicative of opti- 
cally thin synchrotron emission. However, since the map pre- 
sented in Fig. 2 recovers ^ 15 per cent of this integrated flux 
density, we cannot immediately draw the conclusion that 
the imaged structures are also therefore optically thin syn- 



First e-VLBI observations of Cygnus X-3 5 



chrotron emission. Nevertheless, based on simple physical ar- 
guments and comparison with other relativistic jet sources, 
this is likely, and is consistent with the relatively high levels 
of polarisation observed along the jet. As noted above, knot 
C in particular shows a high level of polarisation, which sup- 
ports the interpretation that it is not the core. The highest 
fractional polarisation, almost 25 per cent, appears to the 
east of knot B. This is suggestive of an interaction site be- 
tween the ejected matter and the surrounding medium. The 
same conjecture can be made about the region of relatively 
high polarisation, 15 per cent, in the north of the complex, 
albeit with lower confidence due to the fact that artifacts 
might appear at the edges of the emission. The coincidence 
of emission at the location of knot C in both epochs may 
suggest that this is the location of some standing shock in 
the flow or a site of repeated jet-ISM interactions. 

To summarize, a possible explanation for the radio emis- 
sion detected in the two experiments is the following: the 
radio emitting region of epoch I is a jet-like structure prob- 
ably ejected during the March activity. Epoch II shows mul- 
tiple radio features of what appears to be a one-sided jet. 
The new ejecta interact with the material ejected during 
the previous outburst and generate the high polarisation on 
the eastern edge of the radio structure and a slight bending 
of the jet in this region towards the s outh-west (note that 
precession of the jet is likely a t work - iMioduszewski et all 
l200ll ; IMiller- Jones et all 120041 , thus complicating the inter- 
pretation). The core - corresponding to the location of the 
X-ray binary, and base of the relativistic jet - is probably 
not detected at either epoch, indicating a temporary sup- 
pression of jet production at these times. It is important to 
bear in mind that the jet is ver y close to the line of sight , 
at an angle of less than ~ 15° (|Mioduszewski et a.1.1 l200ll : 
IMiller- Jones et al"1l2004l ) - this will cause an exaggeration of 
apparent bends in the jet, and an underestimate of true jet 
velocities (see more extensive discussions in Mioduszewski 
et al. 2001 and Miller- Jones et al. 2004). 



4 CONCLUSIONS 

We have presented results on the X-ray binary Cygnus X-3 
obtained in the frame-work of the first open call for Euro- 
pean e-VLBI observations. The two successful experiments 
captured Cygnus X-3 initially in an almost quiescent state 
after the outburst in 2006 March, and subsequently in an 
active state a few days after the major flare in 2006 May. 
By the second epoch, a bright, curved, relativistic jet is ob- 
served. The total intensity and polarisation radio maps pro- 
vide evidence for interaction of the jet with the surrounding 
medium and, if our interpretation is correct, with the mate- 
rial ejected in the previous active state. In the first epoch, 
two weeks prior to the onset of the major flare of 2006 May, 
the core is in a deeply radio- quenched state, with a 3a up- 
per limit to the flux density of 0.4 mjy. The use of e-VLBI 
enabled very quick access to the data, practically within one 
day from the end of the runs, compared to sometimes weeks 
in the 'traditional' VLBI experiments. This facility for rapid 
response and analysis of VLBI data is of paramount impor- 
tance for the study of transient sources, enabling decision 
making with respect to potential follow-up observations (not 
only at radio wavelengths) to be much more effective. 
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